Dispersive interactions between atoms and non planar surfaces 
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We calculate the dispersive force between a ground state atom and a non planar surface. We 
present explicit results for a corrugated surface, derived from the scattering approach at first order 
in the corrugation amplitude. A variety of analytical results are derived in different limiting cases, 
including the van der Waals and Casimir-Polder regimes. We compute numerically the exact first- 
order dispersive potential for arbitrary separation distances and corrugation wavelengths, for a 
Rubidium atom on top of a silicon or gold corrugated surface. We discuss in detail the inadequacy of 
the proximity force approximation, and present a simple but adequate approximation for computing 
the potential. 
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I. INTRODUCTION 

Quantum and thermal fluctuations of the electromag- 
netic field in the presence of material boundaries gener- 
ate fluctuating spatial gradients of the field intensity and 
fluctuating induced electric dipoles in ground-state atoms 
close to those surfaces, resulting in an atom-surface inter- 
action via the optical dipole force. These dispersive forces 
have considerable importance in fundamental physics, in- 
cluding deflection of atomic beams close to surfaces m, 
classical and quantum reflection of cold atoms [1, 13| 
and Bose-Einstein condensates (BECs) [H from surfaces, 
dipole oscillations of BECs above dielectric surfaces [l,0|j 
and in possible future applications of single-atom manip- 
ulation in atom chips [1, Q . At the limit of large sepa- 
ration distances (Casimir-Polder limit the interac- 
tions show universal features since they depend only on 
the zero- frequency atom and surface optical responses. 

Dispersive forces can be tailored in different ways, ei- 
ther by eng ineering the optical properties of the surfaces 
pdl . [T^ . Il3j by biasing the temperatures of the surfaces 
and their thermal environment 

0, Q Ell, or by suit- 
ably changing their geometry. Given the complexity of 
these quantum forces, which arise from fluctuations at 
all frequency and length scales, calculations beyond pla- 
nar geometries are exceedingly involved. Until recently, 
approximate methods have been used to deal with non- 
planar setupSjincluding the proximity force approxima- 
tion (PEA) [IJ] and the pairwise summation approach 
(PWS) [131 ■ However, these approximate methods dras- 
tically fail when the surface is not sufficiently smooth 
on the length-scale of the atom-surface separation, since 
they do not properl y ta ke into account the non-additivity 
of dispersion forces [Igl • These large deviations from PEA 
or PWS could be probed by using present-day technol- 
ogy, for example with a Bose-Einstein condensate above 
a microfabricated corrugated surface [l^ . Similar devia- 



tions have been recently measured for the Casimir force 
between a metallic sphere and a nanostructured silicon 
rectangular grating [2O, [2l| . 

A number of methods are available to compute the 
dispersive atom-surface forces (see for example [13, [H, 
m, m, H^)- These methods allow the exact computa- 
tion of the force for very simple geometries, such as an 
atom above a single (possibly multilayercd) plane, sphere 
or cylinder, or an atom inside a plane Eabry-Perot cav- 
ity. Toy models with scalar fields with ideal (Dirichlet) 
boundary conditions have also been considered to com- 
pute the "scalar Casimir-Polder" force for a small sphere 
above a uni-axial corrugated surface p6j . 

Similar problems have been discussed for the Casimir 
interaction between two macroscopic bodies (see for ex- 
ample [23, m, [2^ ) , and the scattering approach has been 
shown to be a very efficient tool for analyzing the case 
of non-trivial geometries [sOl- Eor parallel plates, the 
Casimir force may be written in terms of the reflection 
coefficients seen from the region in-between the plates 
[Sll . [3^ . A similar expression also holds for non-planar 
surfaces, with the specular reflection coefficients replaced 
by more general reflection operators that describe non- 
specular diffraction by the surfaces. This allows for the 
computation of the Casimir force between rough [3^ or 
corrugated [U surfaces. Similar methods have been em- 
ployed to compute the force between a plane and a sphere 
111 m, [13, [33 and between two spheres ^,\^. 

In this paper we develop the scattering approach for 
studying the dispersive interaction between ground-state 
atoms and arbitrarily shaped surfaces characterized by 
frequency-dependent reflection operators. This approach 
provides an exact analytical expression for the two-body 
interaction energy, which can be written in terms of the 
scattering matrices defined for each individual, one-body, 
scatterer. Although this is a simpler problem than solv- 
ing for the many-body Green function, it is still very dif- 



2 



ficult to find the exact scattering matrices for the prob- 
lem of a electromagnetic field impinging on a given non- 
planar geometry. In order to illustrate the method, we 
write the two-body interaction energy in a perturbative 
expansion in terms of the deviation of the surface profile 
from the planar case [l^. Apart from this perturba- 
tive approximation, the approach is valid for any value 
of the corrugation wavelength and includes the full elec- 
tromagnetic fiuctuations and the real material properties 
of the boundaries. We compute the perturbative dis- 
persive interaction potential to first order in the surface 
profile both for ideal and real materials, and discuss the 
prospects of measuring non-trivial (beyond PFA) geo- 
metrical effects with a uni-axial corrugated surface inter- 
acting with a Bose-Einstein condensate used as a vacuum 
field sensor. 

The paper is organized in the following way. In Sec. 2, 
we develop the scattering formalism for the atom-surface 
interaction. The resulting general formula is then applied 
to calculate the potential to first order in the surface 
corrugation amplitude in Sec. 3. We compare our results 
with PFA in Sec. 4 and discuss experimental prospects 
in Sec. 5. Concluding remarks are presented in Sec. 6. 



II. ATOM-SURFACE INTERACTION FOR 
NON-PLANAR GEOMETRIES 

We consider a ground-state atom at position = 
(tatZa) on top of a non-planar surface, as shown by 
Fig. 1. The surface profile is defined by the func- 
tion h{x,y) representing the local height with respect 
to the plane z = 0. The atom is in free space, so that 
ZA > h{rA)- The fluctuating electromagnetic field prop- 
agating from the surface towards the atom is written 
in the mixed Fourier representation (frequency w, two- 
dimensional wave vector k) 

ET(k,z,c.)= ^i(k,c.)e+(k,u;)e^'=-^ (1) 

p=TE,TM 

where p is the polarization index and the longitudinal 
wavevector (sgn denotes the sign function; k = |k|) 

kz = sgn((jj)-\/w^/c2 — k^. 

Similar expressions in terms of unit vectors e-^g and e-^j^j 
hold for the field E^- (k, z, oj) propagating towards the sur- 
face, with the replacement k^ —^ —k^. The complete 
wave vector is denoted as K='= = k± k^z, and the unit 
vectors given by (the frequency dependence will be gen- 
erally omitted) 

4E(k) = e^E(k) = z X k, 
4m (k) = 4E(k) X K±, 

correspond to transverse electric (TE) and magnetic 
(TM) polarizations, respectively. 




X 



FIG. 1: (color online) Ground-state atom on top of a non- 
planar surface (profile function h{x,y)). 

The refiection by the non-planar surface modifies the 
wave vector k and the polarization p = TE, TM, while 
conserving the frequency lu since the surface is assumed 
to be stationary 

EliKu;)^ I ^^(k,p|7^5(c.)|k',p')i^,^.(k',u;). 

p' 

(2) 

The refiection operator TZs{uj) depends on the surface 
profile function h{x,y). Explicit results for its matrix 
elements are presented in Sec. 3 to first order in h{x,y). 

Let us assume for a moment that the z-axis is taken 
along the atom position (ja = 0). At zero temperature, 
the Casimir atom-surface interaction energy is then given 
by the zero-temperature scattering formula as an integral 
over the positive imaginary frequency axis {cu i^) [30| 

C/(Ra) = h / TT^Tr log (1 - 7^5 e-'C^^ TZa e''^'^) , 

/C = diag(K) , K = ^^2/^2 + fc2 (3) 

with TZa and TZs representing the refiection operators 
for the atom and the surface, respectively. /C is a di- 
agonal operator in the basis of plane waves |k,p) with 
eigenvalues k, so that e^'^^-* represents the propagation 
between the two scatterers. Alternatively, it may also be 
interpreted as the displacement operator (soj : here TZa 
is computed for a coordinate axis centered at the atom 
ccnter-of-mass, whereas e~'^^^ TZa e~'^^^ corresponds to 
the 'laboratory' axis. A similar expression holds at fi- 
nite temperature, replacing the integral over imaginary 
frequencies ^ by a sum over Matsubara frequencies. 

It is actually simple to calculate the atomic refiection 
operator for an arbitrary position (tajZa) when taking 
the plane wave basis. The spherically-symmetric ground- 
state atom is described as an induced electric dipole, with 
a dipole moment 

d{Lu) ^ a{uj)E{'RA,i^)- (4) 

Let us emphasize at this point that a{uj) is the dynamic 
polarizability defined according to SI unit system [4]| . 
The dipole field components in the region z < za are 
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written in the representation defined by ([T]) with [42 



e (k) • d{uj)e' 



(5) 



When replacing ^ into (O, we find two separate contri- 
butions associated to the field Fourier components prop- 
agating upwards or downwards. To calculate the atom 
reflection operator, we take the upward components 



(6) 



This result can be cast into a form similar to ^ , with the 
atomic reflection operator given by {lo and —* in) 



(k,p|7^^|k^p') = - 



2k, eoc^ 



Xg— »(k-k')TAg-(K + K')zA^ 

The dependence on ta and za in this equation is ex- 
actly the one expected from the matrix elements of the 
displacement operator in the plane wave basis. The dis- 
placement along the z-axis was already explicitly taken 
into account in the scattering formula ([3]), so we must 
take ZA ~ when replacing ([7]) into ([3]). 

As we assume that the atom-surface separation dis- 
tance is much larger than the atomic dimensions, we may 
expand the general scattering formula ([3]) to first order 
in a{uj) or, equivalently, to first order in TZa [isj 



;7(R^) = -h I 



27: J (2^)2 y (27r) 



d^k' 



p,p' 



Xp\TZs\k',p'){k',p'\nA\k,p) 



(8) 



The ZA dependence is already explicit in ^ because TZa 
in ^ was defined with respect to reference axis such that 
ZA = 0. Note that the resulting za dependence coincides 
with ([7]) , which was derived for an arbitrary position with 
respect to the origin. From ([7]) and we find 



U{Ra) = — / < / JTTT^ 

cqc^ Ja 2n J {2ny 

gi(k— k')-rAg-(K+K')-^A 

2^' 

x^(k,p|7^5|k',p')e+(k)•v(k'). 



d^k' 
?J (27r)2 



(9) 



Note that this general formula holds for magneto- 
dielectric media, including the anisotropic case hj. The 
dispersive potential for an atom with magnetic polariz- 
ability [3| can also be derived from ^ along the lines 
presented above. The corresponding atomic reflection 
operator is calculated in Appendix A. 



As a first application of ([9]), we briefly consider the case 
of a planar surface at z = (corresponding to a profile 
function y) = 0) made of some isotropic material. In 
this case, TZs is diagonal, its matrix elements being 

(k,p|7^s|k^p') = (27r)25(2)(k-k')<5p,p'rP(fc,zO, (10) 

where rP{k,i^) arc the specular reflection coefficients for 
the plane surface. For instance, for a homogeneous non- 
magnetic bulk medium, they are given by the Fresnel 
formulas (e(iC) = electric permittivity) 



r^"(k,ze) ^ , 

K + Kt 



„TM 



(11) 



Eq. ([9]) then yields the known interaction energy between 
a ground-state atom and a plane surface valid for arbi- 
trary separation distances [H, HI, [H, [H, [H] 



eoc^ 



„TE 



27r 



(k,zO 



1-1- 



2c'k 



_fk_ 

(2^ 

,2 7.2 



-2k,za 



2k 



(12) 



„TM 



{k,^0 



This expression can be simplified under the assumption of 
small or large distances za (with respect to some typical 
atomic transition wavelength \a)- For za ^ Xa and 
^ Aa, we get the van der Waals and Casimir- Polder 
limits which correspond respectively to Eqs. (24) and 
(22) in Ref. [H. 

Eq. ([9]) represents a general result for the dispersive 
interaction between an atom and a body scatterer. The 
difficulty, of course, lies in the explicit derivation of the 
matrix elements of the corresponding operator TZs- In 
the next section, we apply this result to derive the lat- 
eral dispersive force for an atom on top of a corrugated 
surface up to first order in the corrugation profile h[v). 
A second interesting application, left for a future publi- 
cation, would be to compute the roughness correction to 
the normal dispersive force - in this case it is necessary 
to compute TZs up to second order in h{v) (see [U 



III. PERTURBATIVE INTERACTION 
POTENTIAL WITHIN THE SCATTERING 
APPROACH 

In order to compute the exact atom-surface interaction 
potential it is necessary to find the reflection operators of 
the surface, which is a highly non-trivial problem. For the 
sake of illustration of the scattering formalism applied to 
the atom-surface problem, we now compute those reflec- 
tion operators in a perturbative expansion in terms of the 
deviations /i(r) of the surface profile with respect to the 
planar configuration. We will calculate the atom-surface 
interaction energy to first order in h. This term gives rise 
both to a correction to the normal force and to a lateral 
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force on the atom which docs not exist in the case of a 
planar surface. We start by expanding the reflection op- 
erators in powers of ft. as 7^ = 7^(") + 7^(l' + 0{h'^). We 
model the optical response of the homogeneous isotropic 
material in terms of its frequency-dependent electric per- 
mittivity e(i^). The zeroth-order term TZ^^'^ is given by 
()10p . whereas the first-order reflection matrix elements 
can be written as 



(k,p|7^W|k^p')-i?i'i(k,k')i^(k-k' 



where H(k.) is the Fourier transform of the surface profile 
h{r). Therefore, the first-order atom-surface interaction 
potential is 



where g^k, za) is the response function 



(14) 



ff(k, Za) 



Ctk' k" 



eoc2 JO 



^-(k +K.")ZA 



ZTT 



dPk' 



r^k'.k'-k 



2k" 



(15) 



^e+(k')-v(k")<j4k',k"). 



The first-order reflection matrix was calculated i n 13311 
by following the approach presented in [i^ (see also [47l|). 
It is a non-diagonal matrix with non-specular reflection 
coefficients given by 



(k, k') = uj,p, {k, k'- iS) A« (k, k'; li), 



(16) 



where we have defined 



tP(k, z^) are the transmission Fresnel coefficients for the 
planar interface defined by 



(13) The matrix A^^^ is expressed as A^-^^ = A^^ — aL^' with 



Ai^)(k,k';zO = (A^t±«:)Br^ 



C S 

S C±l3p[ 



B' 



Bt = diag(l, cKt/ \/ e(jC)0: (i^k/n, f3t = k/Kt, 
C = cos{(j) -(/)'), S' = sin(0-0')- (19) 



We have denoted with (p (0') the angle between k (k') and 
an arbitrarily chosen axis on the surface plane. Starting 
from these definitions and using the scalar product of 
polarization vectors 



e^E(k) • eTE(k') 

^TE(k) • ^TmO^ ) 
^TM(k) ■ ^TeO^ ) 



ck'S 
ckS 



(k).eTM(k') = 



- — {kk' + Kn'C), 



(20) 



Uppi (k, k ; i^) 



rP{k,iC)tP'ik\iO 



(17) we obtain for Ck'^k" the following general expression 



flk'.k" 



-{k' + k.")za ^ 



'wte.te 



wtm.te 



-ute,tm 



[k'k" + K'K"C][t{iC)k'k" 
{eic^){e/c:'~{n'YWO + l]) 



-utm,tm 



(21) 



By inspection of Eq. ((2T|) . one shows that the response 
function g(k, za) depends only on the modulus k = |k| of 
the corrugation wave vector, as expected: the informa- 
tion about the direction associated to the surface profile 
is contained only in H(k), whereas the response function 
derived within the perturbative approach is based on the 
planar geometry and its underlying symmetry. 

The results presented above allow for the numeri- 



cal calculation of the first-order potential for arbitrary 
ground-state atoms and material media by plugging the 
corresponding polarizability and permittivity functions 
into Eqs. ^5]) and ((2T|) . As an illustration, we con- 
sider a ground-state rubidium atom above a gold or a 
silicon surface. The dynamical atomic polarizability is 
obtained from [11], and optical data for gold and silicon 
are obtained from [i^ . The values along the imaginary 
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frequency axis are then calculated with the help of dis- 
persion relations [13, HH ■ 

In Fig. 2, we plot the response function after normal- 
izing it by 



3ficQ;(0) 



(22) 



that is also the result derived by Casimir and Polder 
for the attractive force on an atom on top of a perfectly- 
reflecting plate at large separation distances. The hor- 
izontal axis represents the atom-surface separation za- 
We take two different values of k, corresponding to cor- 
rugation wavelengths A = 27r/fc = 10 /xm and 100 nm. 

The global shape of the curves can be discussed in a 
qualitative manner, which will be made more quantita- 
tive later on. At low values of kzA, PFA is expected to 
hold so that the response is approximately independent 
of k. The ratio g{k^ za) / Fqp thus grows linearly with z, 
showing the well-known power-law change when sweeping 
from the unretarded van der Waals to the Casimir-Polder 
regime. As za approaches 1/fc, PFA becomes gradu- 
ally worse, and then a nontrivial geometry effect reduces 
g{k^ Za) exponentially for za^ ^/k. In the next section, 
we show that geometry and real material effects can be 
disentangled in most cases of interest, thus providing a 
simple manner to explain how they jointly produce the 
results shown in Fig. 2. 
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FIG. 2: (color online) Dimensionless ratio g{k, za) / Fqp{za) 
as a function of separation distance za (for a ground-state Rb 
atom). For the corrugation wavelength A = 2-R/k = 10 /xm, 
we show the values for gold (black) and silicon (blue) surfaces. 
We also show the results for A = 100 nm with gold (dotted 
red) and silicon (dashed green) surfaces. 



IV. DISENTANGLING GEOMETRY AND REAL 
MATERIAL EFFECTS 

As discussed in the Introduction, the most commonly 
used approximation methods to compute atom-surface 



interactions are the proximity force approximation (PFA) 
and pairwise summation approach (PWS). In [l8[ we 
compared these approximate methods with our exact 
scattering approach to first-order in the surface profile 
for the case of perfect reflectors, and showed that PFA 
overestimates and PWS underestimates the lateral atom- 
surface force. In this section we want to expand upon 
these considerations and consider in further detail the 
corrections to the PFA approximation beyond the case 
of perfect reflection. This will allow us to disentangle the 
effects of geometry and real materials, thus leading to a 
very simple rule for the evaluation of g{k, za)- 

In the atom-surface context, the proximity force ap- 
proximation amounts to compute the atom-surface po- 
tential U for a given geometry from the potential for the 
planar geometry taken at the local atom-surface dis- 
tance 

U ~ U^^Hm - h{rA)) =^ U^°^ - h{rA)d,U^'^K 

The PFA holds when the surface is very smooth in the 
scale of the separation distance, that is in the limit kzA — *■ 
0. From our expression for the first-order potential p4|) 
one can indeed prove that for any material the response 
function satisfies the "proximity force theorem" 



g{k^Q,ZA)^~dM^°HzA). 



(23) 



This follows from the fact that the response function 
at zero transverse momentum is given by the specular 
limit of the non-specular coefficients flk'.k", which corre- 
sponds to k' = k". In this case we have i?^} p„ (k', k') = 
2K'rP'(k';iOV,P"' and Ok'.k' - -"^"^-^ ^ 

ep,(k')r-P'(k';ie)- From (US]), 
sponse function is then 



= e-2--E,'<'(k') 
the zero-momentum re- 



g(k = o,zA) = ^ r^ec^m 

eoc^ Jo 27r 

X / ^e-2--^e+(k).e^-(k)r^'(k;*C). (24) 

Comparing with expression (fT2l) for C/^"^ , we immediately 
prove the proximity force theorem (|23p . 

As a consequence of this discussion, we recover ([25)1 
from our more general result (|14p when replacing g(k, za) 
by g{0, za) - Now the value of g{k, za) differs from p(0, za) 
for any finite value of k and it is worth quantifying devi- 
ations from PFA by introducing the function 



pik,ZA) 



9{k,ZA) 
9(0, za) 



(25) 



Deviations from PFA can first be discussed in the vicinity 
of = where PFA is recovered (p ^ 1 for /c — > 0). It 
is possible to deduce from p5 l) - ([2T|) that the first order 
derivative of p with respect to k is identically zero, for any 
atomic polarizability and surface permittivity. It follows 
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that deviations from PFA appear only at order fc^ in a 
Taylor expansion around fc = 



1 



y a^p(0, Za) 



(26) 



For perfect refleetors, we can also prove that the second- 
order derivative in ([26]) is negative, so that p{k^ za) is 
concave near fc = 0. 

When going far beyond PFA (arbitrary values of kzA), 
we have to compute p{k, za) numerically. The result of 
this calculation is plotted on Fig. 3, with p shown as 
a function of kzA for three different values of za in the 
case of an Rb atom on top of a Si surface. We notice that 
these curves are quite close to one another. In order to 
obtain an estimate of their values, we also show the curve 
corresponding to the Casimir-Polder limit for perfect re- 
flectors, which is derived from (j^S]) and the calculations 
of Appendix B 



p^p\k, za) 



-kZA 



1 4- kzA ■ 



16(fczA)2 , {kzAf 



45 



45 



(27) 

For this limiting case, p depends only on the dimension- 
less Z = kzA- We may thus conclude that this vari- 
able captures most of the geometry correction, since p 
depends very little on za for a given Z. For distances 
larger than 1 /im, p may be well approximated by the CP 
formula Eq. ([27|) . This statement holds for the case of 
silicon, but the results for gold surfaces (not shown) are 
even closer to the CP curve shown in Fig. 3. 

Of course, finite conductivity corrections are very im- 
portant for the evaluation of the force between an atom 
and a plane plate, and the full response function g has 
to be written as 



g(fc, z) = p{k, z) riF i^cp 



(28) 



In this formula, F^^p is the Casimir-Polder force ([22|l . 

rjF = -F'-^V'^cp represents the reduction of the disper- 
sive force due to material properties of atom and surface 
(calculated for a plane geometry), and p describes the 
effect of geometry. 

In Fig. 4, we plot -qp as a function of za for a Rb 
atom on top of a silicon or gold surface. As expected, 
the CP result for perfect reflectors is recovered at large 
distances in the case of gold, whereas for silicon we find 
a reduction of rip ~ 2/3 due to the finite value of the 
zero-frequency permittivity in this case. These results 
coincide with those obtained in [i^. For both gold and 
silicon surfaces, the reduction gets stronger as the sep- 
aration distance is decreased, as expected since shorter 
distances correspond to larger frequencies, for which the 
optical responses of both atom and surface are smaller. 
For very short distances, rjp is linear {rjp sa 5.8z(/im) 
for gold and 5.1 z(/xm) for silicon), in agreement with the 
power law modification expected in the van der Waals 
unretardcd regime. 

Since p is well approximated by the Casimir-Polder 
result for perfect reflectors (|27p . the real material and 
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100 nm 
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FIG. 3: (color online) Geometry correction factor p as a func- 
tion of kzA for a silicon surface and za = 100 nm (light blue) , 
1 /im (green), 10 pm (dashed black). The red line corresponds 
to the CP resuh given by Eq. (f27)l . 



beyond-PFA geometry corrections are approximately un- 
correlated. One may compute the former effect for the 
simple plane geometry {r]p{zA)), the latter for 'perfect' 
atom and surface (p^p^kzA)) and flnally combine the 
two effects with the help of (|28p . The simple analytical 
formula (|27|) is thus of great practical relevance, since 
it allows for an easy evaluation of the response function 
g{k, Za)- For a given corrugation wavelength, as was dis- 
cussed in connection with Fig. 2, the two effects are more 
relevant at different separation ranges: r]p is mainly af- 
fected for small za (producing the linear increase of all 
curves on the left part of Fig. 2), whereas p{k,z) dif- 
fers for large kzA (leading to the exponential decay ap- 
parent on the right part of Fig. 2). When k is small 
enough (as in the example with A = 10 pm on Fig. 2), 
the two regions do not overlap, explaining why the max- 
imum can approach unity. Otherwise (as in the example 
with A = 100 nm on Fig. 2), they do overlap and the 
maximum value of the curve remains well below unity. 

Only for very short separation distances (or very short 
corrugation wavelengths for a given kzA) p{k, za) starts 
to deviate slightly from pQp^(fcz^), as illustrated by the 
light blue curve in Fig. 3 for z = 100 nm, resulting in 
some entanglement between real material and geometry 
corrections. In this case, p is well approximated by the 
van der Waals analytical formula derived in Appendix C. 



V. EXPERIMENTAL PROSPECTS 

There are several possible experimental scenarios 
where non-trivial geometrical effects of atom-surface 
forces could be probed. In particular, Bose-Einstein con- 
densates represent an ideal sensor of quantum vacuum 
forces since they are well controlled and characterized. 
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FIG. 4: (color online) Correction r^p = F^°> /F^°^ of the 
normal dispersive force between a Rb ground state atom and 
a plane gold (dark green) or silicon (red) surface as a function 
of the separation distance za- 

One possible experiment is to use a BEG as a sensitive os- 
cillator, whose center-of-mass frequency is modified when 
the BEG is close to a surface (TJ] . This idea was recently 
used to measure the normal component of the Gasimir- 
Polder force between a Rb condensate oscillating along 
the perpendicular direction to a planar dielectric surface 
01 . Lateral Gasimir-Polder forces could be measured in a 
similar fashion with an elongated BEG oscillating parallel 
to the surface, with it long axis parallel to the corruga- 
tion axis of a uni-axial corrugated surface [H, [l^ . We 
have shown in these references that non-trivial geometri- 
cal quantum vacuum effects, beyond the usual proximity 
force and pair-wise summation approximations, could be 
measured using present-day cold atom technology. We 
shall not further discuss this proposal here. 

Rather, let us comment briefly on another possible sce- 
nario, that although is less sensitive, has the nice feature 
of measuring directly the atom-surface interaction poten- 
tial instead of its second derivative, as in the BEG oscil- 
lator. The idea is to use a quasi one-dimensional BEG 
to locally probe potential variations via imaging the lo- 
cal one-dimensional density. Small electric and magnetic 
fields of BEGs in atom chips have been recently probed 
using this technique [1^. Lateral Gasimir-Polder forces 
could be measured with a cigar-shaped BEG placed par- 
allel to the surface and with its long axis perpendicular 
to the uni-axial corrugation lines. The potential along 
the quasi one-dimensional BEG is related to the one- 
dimensional particle density nid{x) as 

Vho(a;) -I- Vcp(x) = -hujtr\/l + 4ascatnid(a;), (29) 

where Vho{x) is the external harmonic trapping potential, 
Vcp{x) is the Gasimir-Polder potential, wtr is the trans- 
verse (strong) trapping frequency, and Cgcat is the s-wave 
scattering length. In [s^ the optimal potential single 



shot sensitivity AT^ is estimated as AV = jAN/ PqXq, 
where 7 = 2h?asciit/ni- Here m is the mass of Rb atoms, 
AiV is the detection imaging noise, xq is the longitudi- 
nal spatial resolution, and po is the transverse spatial 
resolution. In the experiment [s^, Wtr = 27r x SOOHz, 
AiV « 4 atoms per pixel in a GGD camera, leading to a 
single-shot single-point sensitivity to potential variations 
of Ay w 10-13 - 10-"eV. 

For the sake of estimating the order of magnitude of 
the lateral Gasimir-Polder force let us consider a uni- 
axial corrugated perfectly reflecting surface. We now 
use our perturbative expansion of the potential in pow- 
ers of the corrugation profile. In the retarded regime 
{za » A^), the zeroth-order potential is {/^"^(z^) ~ 
-3hca{0)/32n^eoz\, where a(0)/47reo = 47.3 x lO^^Om^ 
for ^''Rb atoms. For a corrugation amplitude h = lOOnm, 
wavelength A = lO/im and an atom-plate distance za = 
2/im, the correction is U^^'' = 1.13 x IQ-^^eV, which is 
on the border of the experimental sensitivity reported 
in {5^ . The experimental signal of the lateral Gasimir- 
Polder force would consist on a density modulation of 
Id BEG density following the law ni^{x) ~ [^o(a;) + 
Vcp{x)]\ 

Other possible experiments to probe lateral Gasimir- 
Polder forces could involve spin echo techniques for 
atomic beams fiying above corrugated surfaces 14], or a 
two-component, phase-separated BEG "level" [SJ] used 
as a cold-atom analog of an AFM to map the corrugated 
surface potential. 



VI. CONCLUSIONS 

We have developed a scattering approach to the dis- 
persive force between a ground state atom and a mate- 
rial body. We have focused on the case of a corrugated 
surface, and derived explicit results to first order in the 
corrugation amplitude. Exact numerical calculation of 
the dispersive first-order potential has been represented 
in terms of the response function g(k, za), calculated for 
arbitrary values of the separation distance za and cor- 
rugation wavelength A ~ 2-K/k, as long as the corruga- 
tion amplitude remains much smaller than both length 
scales. Different atomic species and materials can be con- 
sidered within our formalism. Here we have illustrated 
our method by taking Rb atoms on top of gold or silicon 
surfaces. 

For separation distances larger than 1 /im, the response 
function gik, z) can be easily obtained from the very sim- 
ple GP result (|27p for the beyond-PFA geometry correc- 
tion factor p, together with the numerical values for the 
real material reduction factor rjp calculated for the plane 
geometry. 

The position dependence of the resulting potential 
leads to a lateral dispersive force. We have discussed a 
variety of possible experiments with Bose-Einstein con- 
densates that would allow for the verification of non- 
trivial beyond-PFA effects with current state-of-the-art 
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techniques. 
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APPENDIX A: SCATTERING MATRIX FOR AN 
ATOMIC MAGNETIC DIPOLE MOMENT 

Eq. ([7]) gives the expression of the atomic reflection 
matrix elements under the assumption that the atom is 
characterized by an electrical polarizabihty a^cu). The 
expression of the matrix elements of the operator TZa 
was found starting from the linear relation ([¥]) between 
the induced atomic electric dipole and the external (in- 
coming) electric field calculated in the position and 
writing down the electric field produced by the induce 
dipole moment d(a'). In particular, we were interested 
in the field going towards the surface (J,) as a function of 
the field coming from it (t). 

A completely analogous calculation can be performed 
assuming that the fluctuating magnetic field H induces 
a magnetic dipole moment m on the atom given by 



(Al) 



where is the magnetic polarizabihty. For the field 

generated by the magnetic dipole in the region z < za 
we find the Fourier component 



_g-ik-rAgifciZyl 



xe-(k) • (K- X m{uj)). (A2) 

Using the following expression for the magnetic field com- 
ing upwards from the surface 

HT(k,z,c.) = — ^i?T(k,c.)(K+ X e+(k,c.))e''=^^ 



and replacing Eq. (|Aip into expression (|A2[) we find 



{k,p\n^\k',p') = - 



m) 



-i(k-k')-rA „-(k;+k')z^ 



2k 

xe;(k,zC). K-x (K'+xe+(k',zO 



(A4) 



As a consequence, if we want to include both contribu- 
tions (electric and magnetic dipole moments), the atomic 
reflection operator is simply given by the sum of the oper- 
ators TZa and TZa , given by ([7]) and (jA4p respectively. It 
is possible to check, for example, that for two atoms the 
scattering formula ([3]) gives back the known result for 
the interatomic potential energy in the Casimir-Poldcr 
regime p4| . 



APPENDIX B: CASIMIR-POLDER REGIME 

Simple closed forms for the first-order response func- 
tion can be derived in some limiting cases. In fact, 
g(fc, za) is given by the integral over the positive imagi- 
nary frequency axis of a function involving three quanti- 
ties that decay to zero: the atomic polarizabihty a{iS,)-, 
the round-trip propagation factor exp[— (k' + k")za] and 
the reflection amplitudes rP{k,i^) (due to the frequency 
dependence of the permittivity function). These func- 
tions decay on different frequency scales: c/Aa, c/ za and 
u)p {a, typical frequency associated to the surface optical 
response, for example the plasma frequency for metallic 
media), respectively. As a consequence, the significant 
range of frequencies giving the main contribution to the 
integral in ([TS]) will be determined by the smallest of 
these frequency scales. One can then expect that the 
integral (|15p can be considerably simplified for some par- 
ticular relations between these frequency scales. 

In this Appendix, we consider the case za 2> Ap, Aa, 
usually referred to as Casimir-Polder regime (note that 
we assume za ^ hc/{kBT) so as to neglect thermal cor- 
rections). Since the smallest frequency scale is c/za, we 
can replace the dynamical polarizabihty a(j^) and the 
permittivity e(«^) by their zero frequency values. For the 
metallic case, we can take e(0) —^ oo and then derive 
remarkably simple expressions from pip : 



1 

ak',k" = 2*^ 



-{k' +k")za 



(Bl) 



This expression could also have been obtained by taking 
the usual perfectly-reflecting boundary conditions on the 
corrugated surface. 

After replacing (jBip into 



gcp(fc, za) = 



ha{0) 



di2 
2n 



■ak',k"(zA,C)) 



(A3) we change the integration variables from (^, fc', 0') 



to (k',k",vj), taking k' {k' , k" , ip) 



(k")2 + P]/(2fccos0') and ^(k',k",0') = c[(k')^ - 
(fc')^(K', k", 0')]^/-^. The Jacobian of the transformation 
is I J| =2ck'k"/{4P(k')2cos2 (/)'-[(«:')'- (k")' + A:2]}1/2. 

We first use the property dcf)' {...) = 2 dcj)' {...). 



9 



Then we write the integral as a sum of two contributions: 



dn / dn" 

J\K'-k\ Jo 



and 



dn' / 

Jy/lJU)^+k^ 



dn" 



with cos(/?m = \{>^'Y — W)'^ + fc^]/2K'fc, and introduce 
further auxihary integration variables u = k' + k" and 
V = k' — k" to obtain the final analytical result 



5gf(fc,ZA)=F,^^P^(zA)e-" (1 + Z 
Z = kzA 



AO) 



45 45 



(B2) 



assume an instantaneous atom-surface interaction (limit 
c oo). 

In order to extract some simple expressions for the 
vdW response function, we have first to choose a specific 
expression for the atomic polarizability a{i^). We will 
use the following model function suitable for a multilevel 
isotropic atom having transition frequencies ujno from the 
71-th excited state to the ground state and electric dipole 
matrix elements dno on these transitions 



2 \ " ^nod-nO 



(Cl) 



Then we have to specify the transmission and reflection 
amplitude functions. Here we take Fresnel formulas with 
two different models for the permittivity e(i^). 

Plasma metals. We first take the plasma model for 
metallic materials: 



APPENDIX C: VAN DER WAALS REGIME 

In this appendix, we assume that za ^ Xa- In this van 
der Waals regime, one can neglect retardation effects and 

I 



= 1 + f 



In this case we obtain from (fT5t-(l2T 



(C2) 



;/vdw 



{k,ZA) = - ^ 



192V27reoz3 (a;2 + 2V2xn + 2) 



6V2Z{xn + V2)Ko{Z) + {V2{Z^ + 12)xn+Z^ + 2A)K,{Z)\ , (C3) 



r 



where x„ = ujp/ujno a-nd Z ~ kzA- If we let all the Xn go This result is, as expected, proportional to the plasma 
to infinity, we get the result for a perfect conductor [l9[ frequency ujp. 



kd^ 



? 1927reo4 
\6ZKo{Z) + + 12)A"i(Z) 



(C4) 



In the opposite limit, with ujp <C lo^q for all n, we derive 
from (|C3p the plasmon van der Waals limit 



„plas (-7, ^ , \ _ \^ kd^^Xn 



(C5) 



X [l2ZAo(Z) + + 24)ifi(Z) 



Semiconductors: In order to treat the case of semicon- 
ductors, we assume that the material can be described 
by the Drude-Lorentz model function 



e(*0 = 1 + 



In this case we get 



^I)l(<0)-1) 



■'DL 



e 



(C6) 



gvdwik, ZA 



ikd'^QXn 



UAoZKoiZ) + Ai{Z)Ki{Z] 



„ 384^eo4(7 + 2)5((7 + 2K2 -2)2 
Ao = {j + 2)(v/^x-„ - V2)((7 + 2)xl - 2) 

Ai = 2(7 + 2)i ((Z2 + 12)7 + 24)a;f, - 3%/2(7 + 2){{Z^ + 8)7 + 16)xl - 48(7 + 2)ixn + 2V2{{Z'^ + 24)7 + 48) 

(C7) 
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where 7 = e(0) — 1 and .t„ = cj^i^/ajno- 
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